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OPTIMAL SCALAR PRODUCTS IN THE 
STANDARD LINEAR VISCOELASTIC MODEL* * 

M. PELLICER+ AND J. SOLA-MORALES^ 


Abstract. We study the third order in time linear dissipative wave equation known as the 
Standard Linear Viscoelastic Model, that appears also as the linearization of the so-called 
Moore-Gibson-Thompson equation in Nonlinear Acoustics. We complete the description in 
[5] of the spectrum of the generator of the corresponding group of operators and show that, 
apart from some exceptional values of the parameters, this generator can be made to be a 
normal operator with a new scalar product, with a complete set of orthogonal eigenfunctions. 

Using this property we also obtain sharper decay estimates for the solutions as t —>■ oo, both 
when the operator is normal or not. 

1. Introduction and statement of results 

In this paper we study the third order in time dissipative abstract wave equation 
(1) {u + aut)tt + L{u + I3ut) = 0, with a, /5 > 0 

where L is a self-adjoint, strictly positive operator in a Hilbert space H with compact resol¬ 
vent. Recall that in this situation the eigenvalues fin of L are positive, increasing, tending to 
oo, semi-simple and the corresponding eigenfunctions are an orthonormal family. A typical 
case is when L = —a^A and A is Laplace’s operator with Dirichlet boundary conditions in a 
bounded domain with a regular boundary, and in that case H = 

The dissipative case corresponds io a/(3 < 1. It is known that in this case all solutions 
tend exponentially to zero (see 0. i)- For a/(3 > 1 it can be shown, by looking for the 
appropriate eigenvalues and eigenfunctions or by seeing that the energy is non-decreasing 
(see formula ([I]) below, or 0 . 0 ). that there exist solutions that do not tend to zero when 
t ^ oo. 

When L = —a^A, this equation is called the Standard Linear Viscoelastic model, and 
represents the linear deformations of a viscoelastic solid with an approach that is considered 
to be more realistic than the usual Kelvin-Voigt model. See [5] and [1] and the references 
therein for a discussion of this model. The same equation appears as a linearization for a 
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classical model in nonlinear acoustics, where it is known under the name of Moore-Gibson- 
Thompson equation. See [7] and [H] and their references. Other recent studies on this equation 
can be seen in [2] and [8], for example. The paper [3] deals with the case a//3 > 1. 

The second-order in time strongly damped wave equation, that would correspond to taking 
a = 0 in was studied by the authors in m some years ago, motivated by our previous 
work [To] on a viscoelastic model. In that paper we proved that the inhnitesimal generator of 
the semigroup associated to that equation was self-adjoint in a particular new scalar product, 
that was shown explicitly, provided that the dissipation coefficient was large enough (over¬ 
damping regime). Several optimal decay estimates for the solutions were then obtained by 
using that new scalar product. 

The goal of the present paper is also to show the existence of a suitable scalar product 
associated to the semigroup dehned by ([T|) for some ranges of the parameters. In this scalar 
product the inhnitesimal generator of the semigroup turns out to be a normal operator, hence 
admitting an orthonormal basis of eigenfunctions. The necessary and sufficient condition for 
the existence of this scalar product, and the normality of the operator, is given in Theorem 
m to be stated below. This fact will allow us to obtain the optimal exponential decay rate of 
the solutions. The main point that will have to be checked in the construction of this scalar 
product is that it dehnes a new norm that is equivalent to the natural one. This will be done 
in several functional spaces where the equation ([T]) dehnes a semigroup (in fact, a group). 
The second main result of the present paper is given in Theorem [2J It states the optimal 
decay rate for the solutions, and not only the growth bound of the semigroup, both when the 
scalar product considered in Theorem [1] exists or not. 

Writing U = (u, n, w)^ = (u. Mi, UuY^ we can write ([1]) as the following hrst order evolution 
equation: 


( 2 ) 


— = AU, U e D{A), with AU 

(JjL 


( ^ 

w 



L{u + {3v) 



\ 


The operator A can be dehned in several functional spaces. According to [T] , if L 
with Dirichlet boundary conditions then the hrst possibility is that: 


-a^A 


n = Hliyt) X Hliyt) x 

D{A) = {(m, M, m;), tc G iTg (T2), M-|-/Sm G iT^(T2) n iTg (^)} 


It can be seen that, in this case, {A,D{A)) dehnes a G°-semigroup (in fact a group) and 
that it is dissipative when a < (3. This last part can be seen using the energy functional 
associated to the following scalar product: 

((mi,Mi,m;i), (m2,M2,m;2)) = 

( 4 ) _ _ 

f^(vi + awi){v 2 + aw 2 ) + a? f^V(ui + avi)V(u 2 + av 2 ) + a^a(/3 - a) VviVv 2 
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Observe that if {u,v,w) G H this energy is well defined. From [T] or [5] we know that if 
{u,v,w) G D{A) then 


dE{t) 

dt 


= -a\p-a) / |Vnp 


which exhibits the dissipativeness of the operator when a < (3. 

But this is not the only possible functional framework. According to [9], some possible 
functional settings are: 


-Hi = F)(LV2) X F)(L1A) X H 
n2 = D{L) X D{L) X F)(LV2) 
Us = d\l) X X H 

X D{L) X 


with the corresponding domains for A. Observe that, actually, 'Hi = H defined in ([3]) if 
L = —a? A with Dirichlet boundary conditions. 

The normality of the inhnitessimal generator A in a new explicit metric is given in Theorem 
[U Associated to ([T]) we dehne the following numbers mi, m 2 , which, as it will be seen in 
Section [21 are the zeroes of a certain Cardano discriminant, in fact that of the characteristic 
equation ([8]) below: 


( 6 ) 


with 

(7) 


mi = a- 


-Ci-i/C-2 


8/33 


-, m 2 = a- 


-Ci + y/ C2 

8^3 


Cl 


27-18 




C 2 = Cf - 64 



Theorem 1. If fin 7 ^ mi, m 2 for all the eigenvalues fin of L, then in each of the spaces 
Hi given in (Ej) we can define a new eguivalent and explicit scalar product (l ■)Gi where the 
operator A becomes a normal operator. Also, there exists a set {Tj’*, j = 1, 2, 3, n = 1,..., cxo} 
of eigenfunctions of A which is orthonormal in the corresponding new scalar product and 
complete in Hi. 

Conversely, in the cases where one of the eigenvalues of L coincides with one of the two 
numbers mi or m 2 (including the case mi = m 2 ) then the operator A can not be made to be 
normal in any scalar product. 


Remark 1. Associated to each fin eigenvalue of L, there exist three corresponding eigenvalues 
of A, named A^, A 2 , A 3 , the three solutions of the characteristic equation: 

( 8 ) oA^ T A^ T /3/i^A T fin — 0. 

The role of the numbers mi, m 2 will become more clear in Section [2] (see Proposition Hj) when 
we show that: 

( 1 ) if mi < fin < m 2 , the three of A”, A 2 , Ag are real; 

( 2 ) if fin = mi or fin = m 2 , they are also real, but two of them are equal and not 
semi-simple; 
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(3) if = mi = 1712 then A” = A 2 = Ag G M with algebraic multiplicity equal to three; 

(4) otherwise, A” will be real and Ag = Ag G C \ M. 

None of the situations 1, 2 or 3 was considered to be possible in [9], due to a small error in 
the analysis of the characteristic equation. This error becomes more important if we consider 
a variable and near 0 since in this case the situation /i„ = mg or = ^^2 will happen for 
inhnitely many values of a. 

Remark 2. In contrast with the strongly damped case a = 0 studied in [n] there is no hope 
to obtain that the operator is self-adjoint in any new metric, because of the existence of 
nonreal eigenvalues when « > 0. So, the property of being a normal operator is the best we 
can expect. 

Theorem 2. 

i) Suppose that A is a normal operator in the new sealar product G obtained in TheoremUi 
Then, any solution U{t) of ([2]) decays exponentially in the corresponding norm as 

||t/(f)||G<e"—117(0 )||g, fort>t) 

where (Jmax = (^maxiA) < 0 is the supremum respect to n (which sometimes is a maximum) 
of the real parts of the solutions of the characteristic equation ([H]), that is, the real part of 
the sometimes called dominant spectrum of A (see Proposition\^ below for a description 
^f ^max ) ■ 

a) This decay is optimal in the sense that for each u < a^ax there exist solutions U(t) such 
that 

||17(t)||G —)■ 00 as t —)■ 00 . 

Hi) On the other hand, if we are in the situation where, according to Theorem [21 we can 
not have a new scalar product where the operator A is normal, the previous optimal 
exponential decay rate result still holds in another suitable norm. 

Remark 3. Theorem [2] recalls the decay results of [9], but we slightly improve them in the 
sense that what they prove is that amax is the so-called growth bound of the semigroup, that 

is 

inf{c(; G M; |le"^*|| < Vt > 0} 

(in the usual norm). We show that this inhmum is, in fact, a minimum and also that we 
can take = 1 in a suitable equivalent norm. This is true even for the case when 

(^max = —1//9, the essential spectrum of A (see Propositions E] and E]). 

The results just stated will be developed and proved in the following sections. More con¬ 
cretely, in Section |3] we will show the existence and form of the new equivalent scalar product 
for the given range of the parameters and also prove the normality of the operator in this 
case and the non-normality otherwise (proof of Theorem [1]). Also in Section |3l we will prove 
the optimal exponential decay rate of the solutions in all the cases (that is. Theorem [2]). The 
next Section |2] will be devoted to the description of the spectrum of A, completing the results 
of [9]. 
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2. Description of the spectrum of A 

In this section we make an accurate description of the spectrum of A, (t{A). Most of the 
results can be found in [9], but we include them here for a better global comprehension. 
Nevertheless, there are some differences, that come mainly from the cases where there are 
three real eigenvalues associated to the same value of /i„, that can even be algebraically double 
or triple, cases that were skipped in [9], as it has been said in Remark [TJ Our description is 
summarized in the following three propositions. 


Proposition 3 (The essential spectrum, see [9], Theorem 3.2). In the four functional settings 
considered in ([5]) the essential spectrum of the operator A is 


^ess (-^) 



The definition of essential spectrum can be found in [B] or |1]. 


Proposition 4 (Description of the eigenvalues). The operator A has an infinite number of 
isolated eigenvalues all of them with finite algebraic multiplicity. More concretely, for each 
fin, n E N, eigenvalue of L, there exist three corresponding eigenvalues of A, named A", A 2 
and A 3 , the three solutions of the corresponding characteristic eguation ([H]). Moreover, if 
Aj = On + ibn is nonreal, then On and bn satisfy: 

+ 2a,i (—h PfirA + Tn ( — — ^ 

\a J \a 


(9) 


8 aa^ + Sal 



(10) 6^ — — (Scm^ + 2an + (3fin) ■ 

a 

Under the dissipativeness condition < a < (3, these eigenvalues satisfy the following: 

(1) (a) If ^ ^ < 1, then for all n one of the eigenvalues is real and the other two are 

complex conjugated: A” G M and A 2 = A 3 G C \ M. 

(b) If 0 < ^ < the same happens except, maybe, for a finite number of values of 
n. In this case 0 < ^ < 1, the roots of a certain Cardano discriminant mi, m 2 
given in (j6]) (see the proof for details) are real and satisfy 0 < mi < m 2 . Then, 
if mi < fin < m 2 (this can happen only for a finite number of values of n) then 
Xf, A 2 , A 3 G M. The following exceptional case can also happen (double root case): 
if there exist ni or n 2 such that fi^ = mi or fi^^ = m 2 , then Ag^ = Ag^ G 
M or Xl)^ = Ag^ E is an eigenvalue with algebraic multiplicity egual to two. 
Alternatively, if fin ^ [mi, m 2 ] then A” G M and A 2 = A 3 G C \ M, as in (a). 

(c) If ^ = 1, then mi = m 2 = ^ > 0 and the same as in (a) also happens except if 
there exists ni such that fi^ = In this exceptional case, Xf^ = A^^ = Ag^ = 
— ^EMisan eigenvalue with algebraic multiplicity egual to three. 
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(2) (a) If X is a real eigenvalue of A, then 


( 11 ) 


( 12 ) 


(13) 


(14) 


1 , 1 

a p 


If A is nonreal, then 


^ , < Re{X) < 0. 

p J 

(b) If Pn < Rm and Pn, Rm ^ [mi, m 2 ], then Re^X^) > -Re(A™). 

(c) The limits of the sequences of eigenvalues are the following: 

Imi A" = 

n—>-oo p 

which is not an eigenvalue, but the only element in the essential spectrum of A 
(see Proposition m, 


lim Re{X 2 ) = 

n—)-oo 

and Imi^X^) = C )0 with 


1 1 

a (3 


( 15 ) im{x^) = ^^pi;:+o{pi;:). 

(d) Also, one has that limn^oo R^{X 2 ) is lower than lim„_^oo A” (respectively, equal or 
O' 1 

higher), if — is lower than - (respectively, equal or higher). 


The proof of this Proposition |4] is presented below. The proof focusses in the cases not 
considered in [9]. 

In the next proposition we describe the dominant part of the spectrum, that is the part 
with the highest real part. This real part will be named CmaxiA) (or simply amax when there 
is no confusion). 


Proposition 5 (Dominant spectrum). Let pi be the lowest eigenvalue of L. To find the 
dominant spectrum of A one has to solve the cubic characteristic equation ([H]) with pn = pi- 
Then, 

(1) (a) If the three solutions of this equation are real (including the case of multiple solu¬ 

tions), then the dominant spectrum of A will be {—Ip} and (Jmax = o'maxiP = 
-1/(3. _ 

(b) If the solutions have the form Aj G M and A 2 = A 3 G C \ M, then the dominant 
spectrum of A will be either {A 2 , A 3 } or {-1/(3} (which is not an eigenvalue, but 
the only point in the essential spectrum), or both, depending on which has the 
highest real part. This real part is then amax = OTmaxiA). 

(2) We also claim that all the possibilities can occur as it is shown in the next three 
significative cases: 
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(-^) ^max — ^'^('^ 2 ) (^) — ^'^('^ 2 ) 


< 



X 


< 


(c) (Tmax — P (d) amax — Re(A2) — —1//3 

Figure 1. Plots of the eigenvalues of the operator A (crosses) in the complex 
plane (in solid lines, the real and complex axes), showing different possibilities 

for amax{A). In all of them, the dashed line represents Re(A) = —^ 

which is the limit of the real parts of the nonreal eigenvalues, and the point 
marked as a diamond is — which is the limit of the real ones. 

(a) In the case 1/3 < a/(3 < 1 one is in the situation 1(h) above and the dominant 
spectrum is {A 2 , Ag}. 

(b) IfO < a/(3 < 1/3 and pi is large enough one is in the situation 1(b) above but the 
dominant spectrum is { — 1/(3}. 

(c) If 0 < a/(3 < 1/3 and a is sufficiently small, with fixed (3 and pi, one is in the 
situation 1(a) above and so the dominant spectrum is {—1/(3}. 

In Figure\fiwe can see different examples where the previous situations are attained. 

When the dominant spectrum is {—1/(3}, then there will be no oscillations in the dominant 
part of the solutions and these cases could be qualihed as overdamped. 

Now we proceed to prove the previous propositions. Most of these results can be found 
in [9] but we include the proofs here as they are slight but important differences in some of 
them. 
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Proof, [of Proposition |3] First of all, let /i„, (pn be a fixed eigenvalue and eigenfunction of 
{L,D{L)). We can look for solutions of ([1]) of the form u{x,t) = z{t)(j)n{x). In this case, z{t) 
would be a solution of 

CXZ Z jSfinZ “1“ ~ O’ 

This equation has ([8]) as its characteristic equation. Hence, for each /i„ eigenvalue of (L, D{L)) 
there exist three solutions of (jH]), that form the three sequences A”, A2, A3 of eigenvalues of 
(H,iA(H)), as stated. If any of these solutions is nonreal, we can write it as A” = + ibn, 

with bn 7^ 0. Imposing this in ([H]) and considering separately the real and the imaginary parts 
of the equation, it is easy to see that a„, satisfy (E]) and (ITU]) . 

To see part 1, we need to see whether this three solutions are real or not. For that, we will 
simply apply Cardano’s method to ([H]) for a hxed /i„. The hrst thing to do in this method 
is to apply the change of variable ^ = A + ^ to ([8]) normalized such that the highest degree 
coefficient is equal to one. We obtain 

+ g = 0 


with 

(16) 


P = 




q = 


I~^nl3 qi'n 

ScK^ a 


a 3a^’ ^ 27a^ 

Now, it only remains to look at the sign of the Cardano’s discriminant, that we can think 
as a function of /i„: 


(17) 

with 


d{fin) = + 27g^ = ^d{nn) 


d(^flr 



4 

h'n T n- 


According to Cardano’s method: 

i) if d{fin) > 0, the cubic polynomial has one real root and two complex conjugates 

ii) if d{fin) < 0, the cubic polynomial has three different real roots 

hi) if d{fin) = 0, all of the roots are real, with some of them being multiple. 

Observe that, as > 0, the sign of (lT7l) is determined by the sign of the second degree 

polynomial d{fin)- The roots of this quadratic equation are the constants mi, m 2 dehned in 
(Ej) with the constants Ci, C 2 dehned in ([7]). It can be seen that C2 = 0 if and only of £ = 1 
or ^ = 9. So, depending on whether mi,m2 are real or not and positive (to coincide with 
a value of some /i„), we will have these different possibilities for the sign of the Cardano’s 
discriminant mzi: 

(a) Suppose A < ^ < 1. Then, one can see that C 2 < 0, which means that d{^n) will have a 
constant sign, which is positive. By Cardano’s method, this concludes that ([8]) has one 
real root and the other two are complex conjugates for all n G N. 
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(b) Suppose 0 < I < |. Then, one can see that C 2 > 0 and Ci < 0. This means that 

0 < mi < m 2 are two different positive real roots for d{fin)- This allows us to know the 
sign of d(pn): 

i) d{fin) > 0 if /i„ G (0,mi) U (m 2 , 00 ). So, (|H]) has one real root and two complex 
conjugates if n is such that /i„ G (0,mi) U (m 2 , 00 ) (we have an inhnite number of 
them). 

ii) d{fin) < 0 if pn G (mi, m 2 ). So, ([8]) has three different real roots for any possible n 
such that fin G (mi, m 2 ) (if we have any of them, they will only be a hnite number). 

hi) d{fin) = 0 if there exist rii or 77.2 such that fi^ = mi or = '^ 2 - In this case, 
-^ 2 ^ = or A 2 ^ = Ag^ is a real eigenvalue with algebraic multiplicity equal to two 
(and the other one is a simple real eigenvalue). 

Observe that the last two cases were not considered in [S]. 

(c) Finally, suppose that ^ = |- In this case, in ([7]) we have O 2 = 0 and Oi = —216, so 
mi = m 2 = ^ > 0. This means that d{fin) > 0 for all n G M except for, maybe, the case 
in which there exists fi^ = mi = m 2 , for which the Cardano’s discriminant is zero. But 
also observe that in this situation we also have p = q = 0 (see fllGlB . Hence, by Cardano’s 
method, ([8]) has one real root and the other two are complex conjugates for all n G N 
except for, maybe, this certain ui, for which A"^ = A^^ = Ag^ is a triple real root. If this 
triple root exists, then a simple computation allows to see that it is equal to — which 
is the same as — This case was not either considered in [9]. 


This proves part 1 of Proposition 01 Let us now prove part 2. First, to prove ffTTl) . we write 
as in [9] /(A) = /i(A) + / 2 (A) with /i(A) = aA^ + A^ and / 2 (A) = (3fin^ + Pn, and we recall 
that we are interested in the real solutions of /(A) = 0. It is easy to see that for A < —1/a 
one has /i(A) < 0 and / 2 (A) < 0. Also, /i(A) > 0 for A > —l/jS with /i(—1//9) > 0 and 
/ 2 (A) > 0 for A > —1//3. So, all real roots of /(A) = 0 must be in — 1 /a < A < —1/(3, a.s 
claimed. This is the same argument as that of [9], but we note that also holds in the case of 
three real eigenvalues. 

Suppose now that A is a nonreal root of /(A) = 0. Then, for these values of a, (3 and pn the 
sign of the Cardano discriminant of ([ 8 ]) dehned in flTTl) is positive. It is easy to see that the 
Cardano discriminant of (E]) is (i(/i„)/64, with d{fin) dehned in flTT)) . So, both discriminants 
have the same sign for the same values of the parameters, that is, positive in this case. Then, 
(19|) will also have a single real root, that will be precisely the real part of A we want to bound. 
We dehne g{a) = 8 aa^ + 8a^ + 2 a (A + /3fin) 3- Pn {^ — l), and we write (jH]) as g{Re{X)) = 0. 
We see that 5 f(—l/2(l/a — 1/(3)) < 0 and 5 ^( 0 ) > 0, and applying Bolzano’s Theorem we 
conclude that ffT 2 l) holds. 

To prove claim (b) of part 2 we consider now p as a continuous variable in the open set 
(0,mi) U (m 2 , 00 ) or (0, cxo) depending on mg and m 2 being real or not. According to what 
has been said above, in each of these open sets the equation g{a) =0 has a single and simple 
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real solution a = ao(/i). Deriving implicitly in g{ao{fi)) = 0 we obtain 

dao{fi) _ —(2/3ao(/i) + £ — 1 ) 
dn g'{ao{n)) 

Since we know that ao{fi) > —l/2(l/a —1//3), the numerator is negative. Since the coefficient 
of the cubic term of g{a) is positive and g{a) has only one real root, the derivative of g at 
this root must be positive. 

To conclude the proof of 2(b) we have still to prove that ao{fi) decreases when n jumps 
from p = rrii to /i = m 2 - In more strict words, we want to show that 

lim ao(/i) > lim ao(/i). 

/X— 

Recall that these two values of /x are precisely the values for which the cubic equation 
/(A) = 0 has a double real root, and for mi < /i < m 2 the equation /(A) = 0 has three simple 
real roots, that we can order and call Ai(/x) < A 2 (/x) < A 3 (/x). When /i —)■ mf, two of these 
roots collide, and become precisely ao(mi), and the same happens as /x —)■ m^, when two of 
them collide and become ao(m 2 ). We do not know for the moment which of the three roots 
collide in each case, but it is clear that A 2 (/x) will be involved in the two collisions, since a 
collision between Ai(/x) and A 3 (/x) is not possible without involving A 2 (/x). So, we conclude 
that A 2 (/x) —)■ ao(mi) as /x —)■ mj*" and A 2 (/x) —)■ ao(m 2 ) as /x —)■ m^. 

So, our claim will be proved if we show that A 2 (/x) < 0 for mi < /x < m 2 . The central root 
of a cubic equation /(A) = 0 with three real simple roots and a positive coefficient of the 
cubic term is precisely the unique root that satisfies f'{X) < 0. Then, we can derive implicitly 
with respect to /x in the equation /(A 2 (/x)) = 0 and obtain 

dA 2 (/x) _ -/3A2 (/x) - 1 
dfi /'(A 2 (/x)) 

The numerator is positive because of the upper bound in fllll) . and the denominator is negative 
because of what we just said. This finishes the proof of part 2(b). 

Finally, the proof of part 2(c) can be found in [9] and part 2(d) is a straightforward 
computation. 

□ 


Let us now prove Proposition [5l 

Proof, [of Proposition E] To prove part 1(a) we observe that because of Proposition 0] all 
the real eigenvalues of A satisfy A < —1/(3 and they accumulate at —1/(3. To deal with the 
nonreal eigenvalues we observe that, under the hypotheses of 1 (a), mi,m 2 must be real and 
"Xi < /xi < m 2 . This, together with part 1(a) of Proposition 01 implies that the values of 
fin that will give nonreal roots of ([ 8 ]) will satisfy fin > m 2 . Then, following the proof of 
Proposition 0] part 2(b), and with the same notation, ao{fin) < ao(m 2 ). Even if /x = m 2 is 
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not an eigenvalne of L, the number A = ao{m 2 ) will be a real (multiple) root of of (|H]) with 
1712 in the place of /i„, so the bound flTTD holds for this A, ao{m 2 ) < — 1//5 and 1(a) is proved. 

To prove 1(b) we have just to observe that, as we said, all the real eigenvalues satisfy 
A < —1/(3 and accumulate at this point, and, because of Proposition 0] part 2(b), the real 
parts of the nonreal ones are bounded above by i?e(A 2 ) = i?e(A 3 ). 

To prove 2(a) we observe, as we said in Proposition |4] part 2(d), that if a/(3 > 1/3 then 
the vertical line Re{z) = —l/2(l/a — 1/(3), where the nonreal eigenvalues accumulate (from 
its right hand side), lies at the right of the point z = —1/(3 which is larger than all the real 
eigenvalues, so the real eigenvalues or their limit cannot be dominant. 

Let us prove 2(b). Since a/(3 < 1/3, the point .2 = —1/(3 lies at the right of the vertical 
line mentioned above. Since we know, by (ITT)) , that the function ao(p) as p —?■ oo tends to 
— 1/2(1/a — 1/(3), it is clear that if pi is large enough then ao(pn) < —1//9 — £ for all n and 
some £ > 0. 

To prove part 2(c) we look at the expression of d{^i) as in (ITT)) and observe that d(/ii) < 0 
if a is small enough, so we are in the situation 1(a) and the dominant spectrum will be 

{- 1 // 5 }- 

□ 


3. A NEW SCALAR PRODUCT, THE NORMAL PROPERTY AND DECAY OF SOLUTIONS 

Let be the eigenfunctions and eigenvalues of L: Lcji^ = in ascending order 

(0 < /ii < /i 2 < • • • — oo) and with the collection {</„} being orthonormal in H. We consider 
the associated decompositions of the spaces in ([T]) 1-ii = 0^]^ En, for i = 1,... 4, where the 
En are the three-dimensional spaces spanned by {(</>„, 0, 0), (0, 0), (0, 0, 0n)}- Observe that 

the spaces En do not depend on i and that for all i they are orthogonal to each other with 
the natural scalar products of the spaces "Hj. 

These natural scalar products in the cases i = 1 and i = 3 when restricted to the spaces 
En (and expressed in the previously given basis) are defined by the matrices 



/ ^^n 

0 

o \ 

/ 

0 

0 \ 

(18) 

On,l = 0 

l^n 

0 

, 0„,3 = 0 

t^n 

0 


\ 0 

0 

1 / 

\ 0 

0 

1 / 


We will focus only on 'Hi and H 3 since the spaces H 2 and H 4 can be related to the previous 
ones by the natural isometry 

/ 0 0 \ 

0 0 

\ 0 0 L 02 j 
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Observe also that the scalar product dl]) considered in [T], is different, but equivalent, to 
On,i, as it can be written as 


/ 0 1 0 \ ^ / a^a{P - a)/i„ 0 0 \ / 0 1 0 \ 

I 1 a 0 1 I 0 a^l^n 0 1 I 1 a 0 1 . 

\01a/\ 0 0 

The idea is to define the new scalar product in each of the spaces En by a symmetric 
real matrix Gn,i expressed in the basis {{(j)n, 0 , 0 ), ( 0 , (j)n, 0 ), ( 0 , 0 , 0 ^)} in such a way that the 
eigenfunctions of A, that we call T”’*, T 2 ’* and T 3 ’*, become orthonormal once normalized 
with the natural norm given by On,i- We understand that each Tj’* is the eigenfunction that 
corresponds to A", once expressed in the previous basis. Hence, T"’* = (1, A", (A”)^)^/cJ’* 
where the cj’* is the normalizing constant in the usual norm depending on the space of (jS]). 
Then, we define the matrices 


(19) 


• — 


= col(T”’\ 




and 


-1 

n,i • 


When the matrix Cn,i has the previous form it is easy to see that Gn^i is a real, symmetric 
and positive definite matrix. 

The equivalence between the natural and the new norms is based on the following result. 


Lemma 6. For n sufficiently large, all {x,y,z) G (or C^) and all i = 1,... A there exist 
numbers M,m > 0 (independents of n) such that 


( 20 ) 


M\{x,y,z)\\o„. < 


| g „,. < M\\{x,y,z)\\o„^^ 


Proof. According to Proposition |U for sufficiently large n, the operator A restricted to En 
has eigenvalues A” G M and A 2 = + *6^ = Ag G C \ R. By the same proposition, part 2 (c), 

the limits of A”, and are given by ffT3|) . fflTl) and ([15]). Also, one can easily compute 



and 


h'n T 

+ I + ^ hn + Offin)- 

Then, one can compute the elements of the matrices Gn,i and Gn, 3 . This is not a short 
calculation, and the use of an algebraic manipulator can be helpful. The results, up to 
leading orders as /i„ —)■ 00, are 


( 22 ) 
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/ 2 a ^ J ^ a ±§_ n - Ln(li i /i \ \ 


(23) 


GnA — 


+ 0(/i„) + 0(/in) + o(l) 


2a0 

4a^/l^+^^+3a^ , /i \ 

4«/33 + 


2a 

{a+0)^ 

40/32 


+ 0(1) 


and 


^ + ^(1) 


/ 


/ /i2 + 0(/i2) 


(24) — 




3o/3-o2+y32 _|_ nO/ '1 +±M±+/, -L nO/ 'l 

f^n + 0[^n) 2a8 ^^n + 0[fin) 


fin + 0{fir, 


2o/3 

6a3/3-3a3+2o/32+^3 

4a/3^ 


fif^n “1“ 


/3 

6a^/3-3a^+2a/3^+/3^ i ^ 

-^3- + ®ll^ 


\ 


+ 0 ( 1 ) 


3q^+o/3+/32 

2^2 


+ o(l) / 


Observe that all the leading terms of Gn ,3 are positive because a < j3. Let us now hrst 
prove fl20|) for Gn,i- Observe that, intuitively, this result will be true as both norms have the 
same diagonal terms (in asymptotic order) and the other ones are of lower order so they will 
be controlled by the diagonal ones. To prove that in a rigorous way, consider (x, y, z) G 
(or C^) and consider n large enough such that all the terms of fl25]) are positive (that is, large 
enough such that o(/x„) and o(l) do not affect the sign of the coefficients of the leading terms 
in Gn,i)- Also, we are going to use two inequalities. First: 


(25) - i {c>a- + < .6 < 1 

which is true for any a, 6 G M and c > 0. And secondly, if {sn}, {Pn} are positive real sequences 
such that lim„^oo(sn/'rn) = C > 0, it is easy to see that there exist m,M > 0 such that 

(26) mrn < Sn < Mvn 
if n is sufficiently large. 

Now, we start with the lower inequality of (120)) . From fl2^ and using the left hand side 
inequality of fl2^ 


(27) 


+ 


( 0^+/^ 
1 2a 




+ 0 ( 1 ) 


1 a+/3 
20/3 


+ 0(1)) finV"^ + 


( 0^+^ 
20 


2a0 

1 4o2/32-|-y32_)_3^2 

3 4o/33 


finX^ 


1 (o+/3)2 
Cg 4 q/32 


+ 0(1) 


for certain ci,C2,C3 > 0 such that the previous coefficients are positive for n sufficiently 
large. Observe this is possible just taking ci > 0 and with 1/(3 < cf < and 

C2, C3 > 0 and large enough. Observe also that we can choose these constants independently 
of n. 
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The same idea applies to prove the upper inequality of fl2U]) . From fl2^ and using the right 
hand side inequality of fl2^ . 


(28) 


32- 
2a/32 


+ c| + 0(1) 




«+/3 

2a 


{x,y,z)Gn,i {x,y,z)'^ < 

nil'll n '?f2 I f "+/^ I 1 4 o 2 ^ 2 _,_^ 2 _,_ 3 q ,2 ^ (a+/3)2 

Ofij j finV + ^^+ Ofij ) Z 


1 a+l3 
cJ 2a/3 


for certain C4, C5, Ce > 0 such that the previous coefficients are positive for n sufficiently large. 
In this case this is achieved simply taking C4 = C5 = Cg = 1. 

Finally, with this choice of the constants Cj, (I2B]) holds for (1271) and (125]) . Hence, there exist 
mi. Ml > 0 such that 


mi {ynx‘^ + /Xn2/^ + z'^) <\\(x, y, z )|^ < Mi + fXny^ + z^) 

if n is sufficiently large, which proves this lemma for Gn,i- 

The proof of this result for follows the same idea, but with some slight differences that 
we will point out. Again, consider (x, y, z) E (or C^) and n large enough such that all the 
terms of (12T|) are positive (that is, large enough such that o(/i^), o(yn„) and o(l) do not affect 
the sign of the coefficients of the leading terms in Gn, 3 )- We start with the lower inequality 
of (I2D]) . From (I2T|) and using the left hand side inequality of (125]) . 


{x,y,z)Gn,3{x,y,z)'^ > 


Z2 3al3—a^+/3^ 


(29) 


A'n - + oijJ.D'j 

+ o(/i- 

I (5a‘^+a(3+y^ I n(i: i) 1 6a^P-3a^+2al3'^+y^ , A 

+ [ - W - ^[yf^n + 0[yn)) ^^ - + 0[L)JZ 


+ 


oMa0±^ 
2a/3 Fr; 


1 3aP—a'^+fP ,, 

iif 2^ 


20/32 

Z2 6a^/3—3a^+2o/3^+/3^ 
3 4 q/3^ 


X 

2 


y 


for other Ci, 02,03 > 0 such that the previous coefficients are positive for n sufficiently large 
and of order 0(/in) and 0(1), respectively. For this to be true, we will need to choose 

02 = 02(/in)- It suffices to choose oi > 0, independent of n and such that 0i > ’ 

C2 = C 2 y/JjGi with 02 > 0, independent of n and such that 3a2^"^_)_^2 < C2^ < ^, and 03 > 0, 

independent of n and such that 03 > ^ j 

For the upper inequality, from (I2T|) and using the right hand side inequality of (I25|) one 
gets 


(30) 


X, y, z) On, 3 (a;, //, z)'^ < + 4 ^yn + o(/in)) 

I f <^±oi§±^ I 1 3a/3-a^+/3^ , 2 6Q!^/3-3a^+2a/3^+/33 / - 

'X y 2a0 Fn ^ 2a/32 Fn Og 4o/33 

I /3a2+o/3+/32 I I/O I / I 1 6 q2^-3o3+2o/32+/33 x 


X 

,,2 


for other 04, 05(/in), og > 0 such that the previous coefficients are positive for n sufficiently 
large and of the right order. In this case this is achieved simply taking 04 = Og = 1 and 
C5 i/ y-n- 
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So, as in the case of Gn,i, with the previous choice of the new constants q, (j26l) also holds 
for and flSUD . Hence, there exist m3, M3 > 0 such that 



if n is sufficiently large, which proves the present lemma also for Gn,3- 


□ 


Let us now prove Theorem [TJ 

Proof, [of Theorem [T] When 7^ mi, m2, then in each of the H-invariant three-dimensional 
subspaces En defined in the beginning of this Section there are three different eigenvalues of A 
and one can consider the scalar product given by the real symmetric matrices Gn,i defined in 
([T9]). We can then extend the definition of the scalar product to the whole of l-ii = 0[^i En, 
by a block-diagonal procedure, Gi = diag(Gi, G 2 , ■ ■ ■)■ 

With the scalar product so defined, the subspaces En are orthogonal to each other, so the 
whole set of eigenfunctions 



(31) 


becomes orthonormah The operator A diagonalizes in this basis, its adjoint A* is given by 
just its conjugate matrix, and so A and A* commute and H is a normal operator. 

To see that this new scalar product gives a norm that is equivalent to the old natural 
norm, we use the Lemma [HI above in En^+i © T^no+2 © • • • for no large enough and use in 
El® E 2 ® ... Eno that in finite dimensions all norms are equivalent. 

To hnish the reasoning we have still to prove that the family F* of eigenfunctions is complete 
in each "Hj. Suppose that U' = {u',v',w') G "Hi is a nonzero vector that is orthogonal to the 
whole family F*, and we will arrive to a contradiction. If U' is nonzero, then at least one of 
its three components will be a nonzero element of H. Because of that, it will have at least a 
nonzero component in the basis {(pn), suppose for n = n', so it will have a nonzero projection 
in En', namely Un'- Since the projection U' 1—)■ Un' is orthogonal in all of the cases i = 1, 2, 3 
and 4, we see that (f/', T" = {Un', ’*)Gi5 and this cannot be zero for all j = 1, 2, 3 if 

Un' is nonzero. 

To prove the last part of the Theorem let us suppose = mi < m2 (the cases mi < m2 = 
fin 2 or fini = LUi = m2 are similar). In this case the characteristic equation (|HD has one double 
real root and a different simple real root AiL The restriction of A, as it appears 

in ([2]), to the invariant subspace En^ expressed in the basis {(0ni, 0, 0), (0, pm, 0), (0, 0, Pm)} 
will have the form 



It is easy to see that this matrix has A^^ as an eigenvalue of geometric multiplicity one 
but algebraic multiplicity two. This is a property that will hold independently of the scalar 
product considered. And it is well known that this is impossible for normal operators, that 
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have the property that geometric and algebraic multiplicities of eigenvalues always coincide. 
□ 


Let us now proceed with the proof of Theorem [2l 


Proof, [of Theorem [2] The proof of this theorem is the same in all the spaces given in (|5]). 
Hence, our notation will no distinguish among them and we will not include the superindex 
A which distinguishes among the spaces. 

i) For the parameter values that make A a normal operator in the suitable new scalar 
product G given in Theorem [H it has been shown that there exists an orthonormal and 
complete set of eigenfunctions with HT" = j = 1, 2, 3, n = 1,..., cxo. If 

17(0) = J2nj ^ '^nj and, because of the orthonormality of the 

eigenfunctions. 


iic'Wiil = E s E I'? 


' 711 2i^2i<J'maxi 


(t > 0), 


n,J 
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iii) 


We have seen in Proposition [5] that a^ax is either Re(A2) or —1/(3. In the hrst case, 
the solution U{t) = itself has the optimal decay rate. In the second case, if 

o'max = —1/(3, the sequence A” tends to —1/(3 from the left (see Proposition|H parts 2(a) 
and 2(c)) and the corresponding solutions Unit) = have decay rates A”, which 

can be taken as close as we want to —1/(3. 

The idea of the proof of this part is that when there are non-semisimple eigenvalues they 
cannot be dominant. To proceed in this way, among the sequence 0 < < • • • /^n < 

■■■—)■ oo of eigenvalues of L we distinguish the hnite set S of those that coincide either 
with mi or m 2 dehned in (E]) (see Proposition 0]), parts 1(b) and 1(c)) and accordingly 


decompose H = Hq® Hi and L = 


Lq 0 
0 Li 


in such a way that (t{Lq) = a{L) \ S and 


(t(Li) = S {Hi is hnite dimensional). We make the same corresponding decomposition 


in each of the spaces given in ([5|), "Hj = H/ © Hj, and the operator A = 

Observe that Aq is in the situation described in the hrst part of Theorem [1] and Ai is a 
hnite dimensional operator, with all its eigenvalues being real, and some of them being 
multiple. Hence, according to Theorem [T], we can dehne a new scalar product (•, •)go in 
H/ in which we can obtain the optimal decay inequality of part i) above 


Tlo 0 A 

0 All ■ 


©^“*||go < for f > 0. 


On the other hand, as H\ is hnite-dimensional and according to a well-know result of 
Linear Algebra, for each e > 0 we can dehne a new scalar product (•, Aci.e in Tfi such 
that 

||e"^^*||Gle < QArna,.{Al)+e)t ^ > Q. 
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As it is deduced from Proposition |4] part 2(a), amaxiAi) < —l//d- So, we can choose 
e > 0 such that cXmaxi^i) + £ < —'^/l3 < <Tma3;(^o)- Finally, we define the scalar product 
G[ in Hi as the orthogonal extension of Gq and Gi. It is equivalent to the natural scalar 
product of each Hi because it is so when restricted to each of H^ and H}. And the optimal 
decay rate result follows in the G' norm because the dominant part of the spectrum is in 
(t(Ao) and the optimality is true for Go because of part ii) above. 

□ 
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